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Abstract 



Q \ Corrections of order l/mg [Q = b or c) to the Bethe-Salpeter (B-S) equa- 



tion for Aq are analyzed on the assumption that the heavy baryon Aq is 
composed of a heavy quark and a scalar, light diquark. It is found that in 



' addition to the one B-S scalar function in the limit mq — > oo, two more 

scalar functions are needed at the order l/mg. These can be related to 
the B-S scalar function in the leading order. The six form factors for the 
weak transition Af, Ac are expressed in terms of these wave functions and 
the results are consistent with HQET to order l/mg. Assuming the kernel 
for the B-S equation in the limit mg — > oo to consist of a scalar confine- 
ment term and a one-gluon-exchange term we obtain numerical solutions for 
the B-S wave functions, and hence for the A;, — > Ac form factors to order 
l/mg. Predictions are given for the differential and total decay widths for 
Afe — > AJi^, and also for the nonleptonic decay widths for A;, — > Ac plus a 
pseudoscalar or vector meson, with QCD corrections being also included. 

PACS Numbers: ll.lO.St, 12.39.Hg, 14.20.Mr, 14.20.Lq 



I. Introduction 



Heavy flavor physics provides an important area within which to study many 
important physical phenomena in particle physics, such as the structure and interac- 
tions inside heavy hadrons, the heavy hadron decay mechanism, and the plausibility 
of present nonperturbative QCD models. Heavy baryons have been studied much 
less than heavy mesons, both experimentally and theoretically. However, more ex- 
perimental data for heavy baryons is being accumulated |^, ||, |, H, H and we 
expect that the experimental situation for them will continue to improve in the near 
future. On the theoretical side, heavy quark effective theory (HQET) provides a 
systematic way to study physical processes involving heavy hadrons. With the aid 
of HQET heavy hadron physics is simplified when mq ^ ^qcd- In order to get 
the complete physics, HQET is usually combined with some nonperturbative QCD 
models which deal with dynamics inside heavy hadrons. 

As a formally exact equation to describe the hadronic bound state, the B-S 
equation is an effective method to deal with nonperturbative QCD effects. In fact, 
in combination with HQET, the B-S equation has already been applied to the heavy 
meson system ||^, ^ |TU|. The Isgur-Wise function was calculated |T^] and I/tjiq 
corrections were also considered ||^. In previous work [^, |T2|, T^, we established 



the B-S equations in the heavy quark limit {rriQ — > oo) for the heavy baryons 
Aq and loq^ (where = S, S or ^7 and Q = b or c). These were assumed to 
be composed of a heavy quark, Q, and a light scalar and axial-vector diquark, 
respectively. We found that in the limit uiq oo, the B-S equations for these 
heavy baryons are greatly simplified. For example, only one B-S scalar function is 
needed forAg in this limit. By assuming that the B-S equation's kernel consists 
of a scalar confinement term and a one-gluon-exchange term we gave numerical 
solutions for the B-S wave functions in the covariant instantaneous approximation, 
and consequently applied these solutions to calculate the Isgur-Wise functions for 
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the weak transitions A;, —>■ Ac and ^i*^ ■ 

In reality, the heavy quark mass is not infinite. Therefore, in order to give more 
exact phenomenological predictions we have to include 1 / mq corrections, especially 
l/rric corrections. It is the purpose of the present paper to analyze the l/mq cor- 
rections to the B-S equation for Aq and to give some phenomenological predictions 



for its weak decays. As in the previous work |Tl|, |I2|, we will still assume 

that Aq is composed of a heavy quark and a light, scalar diquark. In this picture, 
the three body system is simplified to a two body system. 

In the framework of HQET, the eigenstate of HQET Lagrangian |Aq)hqet has 
0+ light degrees of freedom. This leads to only one Isgur-Wise function ^(cu) {uj 
is the velocity transfer) for A^ — > Ac in the leading order of the l/rriQ expansion 
p5| , p!6| , [1^, |19|, |2y]. When l/mg corrections are included, another form factor in 
HQET and an unknown fiavor-independent parameter which is defined as the mass 



difference w^Aq — in the heavy quark limit are involved ||T9| . This provides some 
relations among the six form factors for Af, Ac to order l/mq. Consequently, if 
one form factor is determined, the other five form factors can be obtained. 

Here we extend our previous work to solve the B-S equation for Aq to order 
l/mq, in combination with the results of HQET. It can be shown that two B- 
S scalar functions are needed at the order l/mq, in addition to the one scalar 
function in the limit mq — > oo. The relationship among these three scalar functions 
can be found. Therefore, our numerical results for the B-S wave function in the 
order mq oo can be applied directly to obtain the l/mq corrections to the form 
factors for the weak transition A^. It can be shown that the relations among 

all the six form factors for A;, Ac in the B-S approach are consistent with those 
from HQET to order l/mq. We also give phenomenological predictions for the 
differential and total decay widths for A;, — > AJu, and for the nonleptonic decay 
widths for A;, Ac plus a pseudoscalar or vector meson. Since the QCD corrections 
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are comparable with the l/mq corrections, we also include QCD corrections in our 
predictions. Furthermore, we discuss the dependence of our results on the various 
input parameters in our model, and present the comparison of our results with those 
of other models. 

The remainder of this paper is organized as follows. In Section II we discuss the 
B-S equation for the heavy quark and light scalar diquark system to order l/mg 
and introduce the two B-S scalar functions appearing at this order. We also discuss 
the constraint on the form of the kernel. In Section III we express the six form 
factors for A;, — Ac in terms of the B-S wave function. The consistency of our 
model with HQET is discussed. We also present numerical solutions for these form 
factors. In Section VI we apply the solutions for the A^ —>■ Ac form factors, with 
QCD corrections being included, to the semileptonic decay A;, AcW, and the 
nonleptonic decays A;, — > Ac plus a pseudoscalar or vector meson. Finally, Section 
VI contains a summary and discussion. 

II. The B-S equation for Ag to l/mg 

Based on the picture that Aq is a bound state of a heavy quark and a light, 
scalar diquark, its B-S wave function is defined as pJ] 



where ipQ^Xi) and (p{x2) are the field operators for the heavy quark Q and the 
light, scalar diquark, respectively, P = m\gV is the total momentum of Aq and 
V is its velocity. Let mq and be the masses of the heavy quark and the light 
diquark in Aq, p be the relative momentum of the two constituents, and define 



X{XI,X2,P) = (0|T^q(xi)v9(x2)|Aq(P)), 



(1) 



A 



1 — 



rriQ+mi:, 



rriQ+mn 



. The B-S wave function in momentum space is defined 



as 
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where X = XiXi + X2X2 is the coordinate of the center of mass and x = xi — X2- 
The momentum of the heavy quark is pi = XiP + p and that of the diquark is 



P2 = — A2-P + p. Xp{p) satisfies the following B-S equation 

Xp{p) = Sp{X^P + p) J ^^K{P,p,q)xp{q)So{-X2P + p), (3) 



where K{P,p,q) is the kernel, which is defined as the sum of all the two particle 
irreducible diagrams with respect to the heavy quark and the light diquark. For 
convenience, in the following we use the variables 



Pi = vp- AamAo , pt=p - {v ■ p)v. 



(4) 



It should be noted that pi and pt are of the order Aqcd- The mass of Aq can be 
written in the following form with respect to the l/mq expansion (from HQET): 

1 



= mg + mz) + H Ei + 0(l/mn), 

niQ ^ 



(5) 



where £"0 and Ei / rriQ are binding energies at the leading and first order in the 1 / mq 
expansion, respectively. rriD, Eq and Ei are independent of mq. 

Since we are considering I/tjiq corrections to the B-S equation, we expand the 
heavy quark propagator Sf{XiP + p) to order l/mq. We find 



Sp — Sqf H SiF, 

mq 



(6) 



where Sqf is the propagator in the limit mq 00 

1 + ^ 



OF 



2{pi + Eo + niD + ie) 



(7) 



and 



S 



IF = I 



;-Ei+p?/2)(i + ^) 



+ 



2{pi + Eo + ruD + ief ' 2ij)i + Eq + nio + ie) ~ 
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It can be shown that the hght diquark propagator to l/mq still keeps its form 
in the limit mq —>■ oo, 

Sd = — T;r.g , . , (9) 



where Wp = ^Jpf + m|,. 

Similarly to Eq.(P), we write Xp{p) K{P,p, q) in the following form (to order 

XP{P) = Xop{p) + —X1P{P), K{P,p,q) = Ko{P,p,q) + —K,{P,p,q), (10) 

mq 

where Xip{p) Ki{P,p, q) arise from 1/mq corrections. As in our previous work, 
we assume the kernel contains a scalar confinement term and a one-gluon-exchange 
term. Hence we have 

-iKo = I®IVi + v^®{p2+p'2YV2, 

-iKi = l0lV3 + -f^0ip2+P2rV,, (11) 

where in Kq appears because of the heavy quark symmetry. 

Substituting Eqs.(|^) and (0) into the B-S equation we have the integral 
equations for Xop{p) Xip{p) 

Xopip) = SoFiXiP + p) I -^^MP,p,q)xopiq)SDi-\2P + p), (12) 



and 



Xipip) = SoFiXiP + p)J-0^,K,iP,p,q)xopiq)SDi-X2P + p) 
+S,f{X,P + p) J -0^KoiP,p,q)xopiq)SD{-X2P+p) 
+SoFiXiP + p) I -^Ko{P,p,q)xipiq)SDi~X2P + P). (13) 
Eq.(^) is what we obtained in the limit mq oo, which together with Eq.(|^) gives 

i'Xopip) = Xop{p), (14) 



since ff = = \ and so fS^p = Sqf- Therefore, Sof{XiP +p)'yfiXop{(l) = Sof{^iP+ 
P)ViiXop{<l) in the first term of Eq.(|T3|). So to order l/mg, the Dirac matrix 7^ from 
the one-gluon-exchange term in Ki{P,p, q) can still be replaced by v^. 
We divide Xipip) i^ito two parts by defining 

Xip{p) = xtp{p) + Xip{p), fXip{p) = ±X^p{p), (15) 

i-e., xtpip) = Uxipip) + i'Xipip)] and Xipip) = ^[Xipip) - fxipi.p)]- After writing 
down all the possible terms for Xop{p) aiid Xipip)-, a^id considering the constraints 
on them, Eqs.(0) and (|T5|), we obtain that 

Xop(p) = 0op(p)maq(^^,s), 
xtp{p) = (I^ip{p)uaq{v,s), 

Xip{p) = hp{p)MAQ{v,s), (16) 

where 4>op{p), 0ip(p) and 4>2p{p) are Lorentz scalar functions. 

Substituting Eq.(|l|) into Eqs.(|l|)(0) and using Eqs.(0)(D(D we have 

1 f d^q 

<Pgp{p) = -7 , . w 2 — un I • A / 7T^^o(^,P,g)0op(g), (17) 

[pi + Eo + mp, + ie)[pf - + le) J (27r)4 



1 f d^ij 

[pi + Eq + mo + ie){pf - + It) J (27r)4 



1 /• d^o 



(pi + + mo + ieM - Wl + «e) 7 (27r)4 

0op(g), 



and 



02P(P) = ^0OP(P)- (19) 

t>op{p) is the B-S scalar function in the leading order of the l/mg expansion. 



which was calculated in [|TT|. From Eq.(|T9D 02p(p) can be given in terms of (popip) 
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The numerical solutions for (f)op{p) and <f)ip{p) can be obtained by discretizing the 
integration region into n pieces (with n sufficiently large). In this way, the integral 
equations become matrix equations and the B-S scalar functions 4>op{p) and (pip{p) 
become n dimensional vectors. Thus 0op(p) is the solution of the eigenvalue equation 
{A — I)(j)Q = 0, where A is an n x n matrix corresponding to the right hand side of 
Eq.(|l^). In order to have a unique solution for the ground state, the rank of {A — I) 
should be n — 1. From Eq.(|l8]), (pip{p) is the solution of {A — I)(f)i = B, where B is 
an n dimensional vector corresponding to the second integral term on the right hand 
side of Eq. (p!8D . In order to have solutions for 0ip(p), the rank of the augmented 
matrix {A — I,B) should be equal to that of {A — I), i.e., B can be expressed as 
linear combination of the n — 1 linearly independent columns in {A — I). This is 
difficult to guarantee if B ^ 0, since the way to divide {A — I) into n columns is 
arbitrary. Therefore, we demand the following condition in order to have solutions 
for 0ip(p) 



yf?! W iP, P, <l) + ^\. Ko (P, p, q) 



ct)op{q) = 0. (20) 



Eq.(^) provides a constraint on the form of the kernel Ki{P,p,q), in which Ei is 
also related Ki{P,p,q). In this way, (pipip) satisfies the same eigenvalue equation 
as (pQpip). Therefore, we have 

0ip(p) = (r(t)Qp{,p), (21) 

where a is a constant of proportionality, with mass dimension, which can be deter- 
mined by Luke's theorem p2| at the zero-recoil point in HQET. We will discuss it 



in the next section. 

Since both (pipip) and (p2p{p) can be related to (/'op(p), we can calculate the 
l/rriQ corrections without explicitly solving the integral equations for (j)ip{p) and 
<p2p{p)- In the previous work [|TT| (f)Qp{p) was solved by assuming that Vi and V2 in 



Eq.(^) arise from linear confinement and one-gluon-exchange terms, respectively. 
In the covariant instantaneous approximation, Vi = yi|p,=g, , i = 1,2, we find 

^1 = 77 , - i27iy6%pt - qt) 



^771 r\^^~^777I „X2 , 021' l^^^^ 



l_67r a(^ff)^gg 

wfiere k and a^'^''^^ are coupling parameters related to scalar confinement and the 
one-gluon-exchange diagram, respectively. They can be related to each other when 
we solve the eigenvalue equation for (pop{p). The parameter /x is introduced to avoid 
the infra-red divergence in numerical calculations, and the limit /i — is taken 
in the end. It should be noted that in V2 we introduced an effective form factor, 

{off)„2 

F{Q'^) = "q2^q-2 , to describe the internal structure of the light diquark [ p3| . 
Defining 0op(pt) = / ^<Pop{p) the B-S equation for (popipt) is Q] 

0op(p*) = ^ / 7^(^i - 2W,vSopiqt), (23) 

in the covariant instantaneous approximation. The numerical results for 0op(^t) can 
be obtained from Eq. (P^D , with the overall normalization constant being fixed by the 



normalization of the Isgur-Wise function at the zero- recoil point [Til]. Furthermore, 
(pop{p) is expressed in terms of (pop{qt)- 

<l^opip) = , , ^ , , \, 2 .... , • ^ / 7^(^i + ^PiV2)4>opiqt). (24) 



{pi + Eo + mD + ie){pf - + ie) J (27r)= 
III. A5 Ac form factors to l/mq 

In this section we will express the six form factors for the A;, Ac weak transition 
in terms of the B-S wave function and show the consistency between our model and 
HQET. 

On the grounds of Lorentz invariance, the matrix element for Af, A^ can be 
expressed as 
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-(Gi(a;)7^ + G2{uj)v, + G,{uj)v'Mu^,{v), (25) 



where is the V — A weak current, v and v' are the velocities of A;, and Ac, 
respectively, and uj = v' ■ v. 

The form factors Fi and (i = 1, 2, 3) are related to each other by the following 
equations, to order l/rriQ, when HQET is applied 



F2 



G2 — —Gi 



rric rrih/ I + uj 
1 A 



— G3 — — Gi" 



rric 1 + uj' 
1 A 



(26) 



mf, 1 + a; ■ 

where A is an unknown parameter which is defined as the mass difference — ttlq 
in the limit niQ 00. 

On the other hand, the transition matrix element of A^ —>■ A^ is related to the 
B-S wave functions of A5 and A^ by the following equation 



{Ac{v')\J,\M{v)) 



(27) 



where P {P') is the momentum of A;, (A^). Xp'{p') is the wave function of the final 
state Ac(f') which can also be expressed in terms of the three B-S scalar functions 
0op(p), 0ip(p) and 02p(p) in Eq.(||) 



Xp{p) = uaq{v, s) \ (popip) + [0ip(p) + 4>2p{p)A] 



(28) 



Substituting Eqs.(|^) and (|^) into Eq.(^) and using the relations in Eq.(^) 
we find the following results by comparing the 7^, 7;_t75, f^(l — 75) and f^(l + 75) 
terms, respectively: 



rric rribJ 1 + UJ 



A 



0^,{<PoP'{k')<Pop{k){kf-W',] 



+ — [0iP'(^') - m + mn)(l)2P'{k')](l)opik)ikf - Wl 
rrir 



+— (-/i + /2 + 2mDF) + — (/i - /2) 
rUc nib 

+—(P,p,{k')[Mk) - {h + mD)ci>2p{mk'i - Wl] 



+0{l/ml), 



(29) 



—I 



(27r) 



nib 



+ —/oP'ik')[Mk) - ik + niD)hpmikf - Wl) \ + 0(l/m^), (30) 



1 



1 



A 



1 + a; 



iGi- 



2(/i - mz)F) 



o(iK 



A 



2/2 



0(l/m^), 



(31) 



(32) 



l + UJ 

where we have defined /i, /2 and F by the following equations, on the grounds of 
Lorentz invariance: 

d^k 



(2vr)^ 



:(t>2Ak')(^op{k){k^-Wl) = F, 



d^k 
(2^ 



ct>2P'{k')cl>op{k)k^^{k^i - wl) = hv^ + hv'^. 



(33) 
(34) 



Eq.(ll) leads to 
/1 + /2 



1 



d^k 



l + ujJ (27r) 



(t^2P'{k')cPopmkl - Wl){v -k + v'-k). (35) 



Eqs.(^9]) and (|30D give the expression for Gi to order l/mg. From Eqs.(pID and 
(|3^ ) we can see that Eq.(|2^) is the same as Eq. (|30|) . Therefore, we can calculate 
Gi to l/ntQ from either of these two equations. This indicates that our model is 
consistent with HQET to order l/mq. 

Substituting Eq.(|55|) into Eq.(^D|) and using Eq.(^) we have 

d^A; 



Gi 



j^,[ci>oP'{k')cl>op{k){k', - Wl) 
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+—^,p'{k')[Mk) - \{h + mD)<Popm{kl - Wl) 

Tfll, Z 

+ (— + —) 7^77^Mp'{k')Mk){kf - W',){v .k + v'-k)]. (36) 
Kmc mb/ 2(1 + u) J 

The first term in Eq.(^) gives the Isgur-Wise function which was calculated in 
our earlier work ||lT|. In order to obtain the l/mq corrections, we have to fix (f)ip{k). 
Fortunately, this can be done by applying Luke's theorem The conservation 
of vector current in the case of equal masses for the initial and final heavy quarks 
leads to 

G,{u = l) = l + 0{l/m%). (37) 
Thus from Eqs.(^), (^) and ( |T7D we have 

a = 0. (38) 

Therefore, (pipik) does not contribute to Gi. 

Now we calculate Gi through Eq. (p6D . Since in the weak transition the diquark 
acts as a spectator, its momentum in the initial and final baryons should be the 
same, p2 = v'l- Then we can show that to order l/mq 

k[v' + k[ = kiv + kt. (39) 

From Eq. (|39|) we can obtain relations between /c[, k[ and /c;, kt straightforwardly: 
k'l = kiuj — kt\/uj'^ — lcos6', 

fc|2 ^ ^2 ^ ^2^^2 _ ly^^^Q ^ ^2^^2 _ 1) _ 2kiktUjVuj^ - Icos^, (40) 

where 6 is defined as the angle between kt and f^. 

Substituting the relation between (pop{p) and 4>op{Pt) [Eq.(p^)] into Eq. (p6D , using 
the B-S equation (P5|), and integrating the ki component by selecting the proper 
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contour we have 



GiH = e(^) + —A,{u) + — AfcH, (41) 
nic nil, 



where 



eM = -/|^^K^*)' (42) 
A(-) = - / ^ ^,), (43) 



, , , , d^kt kt^/oj"^ — lcos6' ^ , , , 



and F{uj, kt) is defined as 



Eo + niD - ujWk - ktVu^ - lcos6' ^ (27r)^ 



[Fi(fcJ - n) - 2{ujWk + ktVU^cose)V2{K - n)] \k,=-w, ■ (45) 



The three dimensional integrations in Eqs. ( ^^ - |4^) can be reduced to one dimen- 



sional integrations by using the following identities: 

d^g^ p{qj) fq^dqt 2p{qf 



(27r)3 [(p, - q,f + ^2]2 J 4^2 (^2 + ^2 + ^2)2 _ 4p2^2 ' 



(46) 



and 



/i_2 ol„ ll„ l^^/l„ I I r5i' \^') 





+ \qt\? + 5' 









(27r)3(p,-g,)2 + 52 

where p{q}) is some arbitrary function of q^. 

In our model we have several parameters, a^s^\ k, Qq, mn, Eq and Ei. The 
parameter Qq can be chosen as 3.2GeV^ from the data for the electromagnetic form 
factor of the proton [^. As discussed in Ref.[^, we let k vary in the region between 



0.02GeV^ and O.lGeV'^. In HQET, the binding energies should satisfy the constraint 
Eq. (!]) . Note that itlo + Eq and Ei are independent of the fiavor of the heavy quark. 
From the B-S equation solutions in the meson case, it has been found that the values 
nil, = 5.02GeV and = 1.58GeV give predictions which are in good agreement 
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with experiments [Q. Since in the b-baryon case the 0(l/m^) corrections are very 
small, we use the following equation to discuss the relations among mD, and Ei, 

mD + Eo + —El = O.Q2GeV, (48) 
ruh 

where we have used m-A^ = 5.64GeV. The parameter rriD cannot be determined, 
although there are suggestions from the analysis of valence structure functions that 
it should be around 0.7GeV for non-strange scalar diquarks Hence we let it vary 
within some reasonable range, 0.65GeV ~ 0.75GeV. In the expansion with respect 
to the heavy quark mass, we roug hly expect {^Ei)/Eo ~ Therefore, Ei 

should be of the order Aqcd-^'o- In our numerical calculations, we let /3(= Ei/Eq) 
change between 0.2 and 1.0. Then for some values of m^) and j3 we can determine 
Eq. Using Eqs. ( p^]4^) and Eq.(p6|) we obtain numerical results for the weak decay 



form factors Fj, Gi {i = 1,2,3) to order l/rriQ. It turns out that the numerical 
results are very insensitive to the value of j3, so we ignore this dependence. We 
also find that the dependence of Fi, Gi {i = 1, 2, 3) on the diquark mass mo is not 
strong. In Fig.l we plot the numerical results for Fi {i = 1,2,3) for k, = 0.02GeV^ 
and K = O.lOGeV^, respectively, with rrif) = 0.7GeV. 

IV. Applications to AJi/ and A^, AcP{V) 

With the numerical results for Fi, Gi {i = 1,2,3) to l/mq obtained in Sec. Ill, 
we can predict the A5 semileptonic and nonleptonic weak decay widths to 

order I/uiq. Since the QCD corrections to these form factors are comparable with 
the l/mq effects, we will include both of them to give phenomenological predictions. 



Neubert |25] has shown that the QCD corrections to the weak decay form factors 



can be written in the following form (up to corrections of the order asA/mg): 

TT TT 

AF, = -^^^v., AG. = -e^^a„ (^ = 2,3), (49) 

TT TT 
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where Vi = Vi(uj) and a,- 



1, 2, 3) are the QCD corrections calculated 



from the next-to-leading order renormalization group improved perturbation theory. 
The scale m is chosen such that higher-order terms (asln(mfc/mc))"' (n > 1) do 
not contribute. Consequently, it is not necessary to apply a renormalization group 
summation as far as only numerical evaluations are concerned. It is shown that m 
can be chosen as 2mi,mc/{mi, + rric) — 2.3GeV. The detailed formulae for Vi and Oj 



can be found in ||2^, which also includes a discussion on the infra-red cutoff employed 



in the calculation of the vertex corrections. As in 25 , we choose this cutoff to be 



200MeV which is a fictitious gluon mass. Furthermore, we use Aqcd = 200MeV in 
our numerical calculations. 
1.4 
1.2 




UJ 



Fig.l The numerical results for Fi {i = 1,2,3) for k = 0.02GeV^ (solid lines) and k = 
O.lOGeV^ (dotted lines), with m/j = O.TGeV. From top to bottom we have Fi, F3, and 
F2, respectively. 



A. Semileptonic decays 



Making use of the general kinematical formulae by Koner and Kramer PB[, we 
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find for tlie differential decay widtli of A;, 
dr 2 



-m 



-A (— + —] [3{r] + 7]-') - 4 - 2cu] + — t;i(cu - l)[3{r] + r]-^) +2-iu] 
' mi)/ 



rrir 



TT 



a. 



+— ai(cu + l)[3(r/ + r/-^) - 2 - Alu] - —{u' - l)[t;2(l + r^) + v^il + 



-1^ 



TT 



77 



+a2(l-r7) + a3(r/-^- 1)]}, 



(50) 



where r] = m\^/miy^ and A = j^^\Vcb\'^ B {K^ ah), witli \Vch\ being tlie Kobayasfii- 
Maskawa matrix element. i?(Ac ^ a6) is the branching ratio for the decay A^ — >■ 



■1 + ^ 



+ 6(0 ) through which Ac is detected, since the structure for such decays is 



aheady well known. It should be noted that in Eq.(^^ ©(ctsA/mg) corrections have 
been ignored and the lepton mass is set to zero. The plot for ^^^^ is shown in 
Fig.2 for rriD = 700MeV, where we also show explicitly the effects of both l/mq 
and QCD corrections. For other values of the results change only a little. 

120 



100 - 



(GeV^) 60 




1 1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.4 1.45 

Fig.2 The numerical results for A~^^ for k = 0.02GeV^ (solid lines) and k = O.lOGeV^ 
(dotted lines), with mo = 0.7GeV. From top to bottom we have the predictions without 
l/niQ and QCD corrections, with l/mq corrections, and with both l/mq and QCD 
corrections, respectively. 
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Table 1: Predictions for the decay rates for A;, AJz/, in units 10^'^s ^B{Ac ah) 



mB(GeV) 


To 




Tl/mQ+QCD 


0.65 


4.77 (7.20) 


4.26 (6.62) 


3.10 (4.76) 


0.70 


5.12 (7.12) 


4.60 (6.56) 


3.34 (4.72) 


0.75 


5.40 (7.02) 


4.89 (6.50) 


3.54 (4.67) 



After integrating uj in eq. (^) we have the total decay width for A;, — > KJ,!). The 
numerical results are shown in Table 1 for m£, = 650MeV, 700MeV, 750MeV and for 
K = 0.02GeV'^ (O.lOGeV^). Fq, Ti/mQ and Fi/mQ+QCD are the decay widths without 
1/mQ and QCD corrections, with l/mg corrections, and with both l/mq and QCD 
corrections, respectively. We have used Vcb = 0.045 in the numerical calculations. 

We can see from Fig. 2 and Table 1 that both l/mq and QCD corrections reduce 
the decay width for A^ —>■ AJiy, and the QCD effects are even bigger. From Table 1 
we can also see that the dependence of our predictions on is not strong. Q 

B. Nonleptonic decays A5 ^ Ac P {V) 

In this subsection we will apply the numerical solutions for the form factors 
Fi, Gi {i = 1,2,3) to the nonleptonic decays A^ AcP{V) {P and V stand for 
pseudoscalar and vector mesons respectively). The Hamiltonian describing such 
decays reads 

H,s = ^V,tV{jr>{aiOi + a202), (51) 

with Oi = {DU){cb) and O2 = {cU){Db), where U and D are the fields for light 
quarks involved in the decay, and (^1^2) = ^?i7/^(l — 15)<12 is understood. The pa- 
rameters Oi and 02 are treated as free parameters since they involve hadronization 

effects. Since A^ decays are energetic, the factorization assumption is applied so 
"'^We note that the results without either l/mq and QCD corrections in Table 1 are bigger than 
those presented in Ref. by about 18%. This is because we employed a cutoff in the numerical 
integrations in Ref. ||ll|] , while the integrations are carried out to infinity in the present work. 
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that one of the currents in the Hamiltonian ( |51| ) is factorized out and generates a 
meson[^, Thus the decay amphtude of the two body nonleptonic decay be- 
comes the product of two matrix elements, one is related to the decay constant of 
the factorized meson (P or V) and the other is the weak transition matrix element 
between Af, and 

M^-(A, ^ A,P{V)) = ^K.l^^^ai(P(\/)|A,(V;.)|0)(A,(POI^1A.(i^)), (52) 

where (0|y4^(V^)|P(V^)) are related to the decay constants of the pseudoscalar meson 
or vector meson by 

{0\A^\P) = ifpq,, 

{m,\V) = fvmve,, (53) 

where is the momentum of the meson emitted from the W-boson and is the 
polarization vector of the emitted vector meson. It is noted that in the two-body 
nonleptonic weak decays A;, AcP{V) there is no contribution from the 02 term 
since such a term corresponds to the transition of A;, to a light baryon instead of Ac. 
On the other hand, the general form for the amplitudes of AcP{V) are 

M(A,^A,P) = tu^^{P'){A + B-f,)uA^{P), 

M{A,^A,V) = nA.(P')e*'^(Ai7^75 + v42P;75 + 5i7^ + 52P>A,(P). 

(54) 

Alternatively, the matrix element for A;, — > A^ can be expressed as the following 
on the ground of Lorentz invariance 

(A,(P')|J^|A,(P)) = UAAn[fl{q'h, + ^f2{q')a,.q'' + f3{q')q, 

-{giiq^hti + i92{q'^)(^^uq'' + 93{q^)qtih5\uA,{P), (55) 

where Qi [i = 1,2,3) are the Lorentz scalars. The relations between Qi and 

17 



Fi, Gi are 



/i 

h 

91 

92 
93 



^2 -F3 



2 \mA^ mA, 
1 ( F2 F3 



2 \mA6 mA,/ 
1 / G-j G-i \ 



G2 G3 



2 \mA, mA, 
1 ( G2 G3 



(56) 



2 \mA, rriKj 

The decay widths and the up-down asymmetries for — > KcPiV) are available 
in Refs.pl 



r(A6 ^ A,p) 

«(A5 A,P) 



Stt 



(mA, + mAj2-m^ 



P\ A\2 



m 



\A\^ 



(mA,-mAj2-m 



-P I Dl2 



151^ 



Ab 



2\P'\Re{A*B) 



where A and B are related to the form factors by 



(57) 



A 
B 



G 

^ycbVuo(^ifp[{mA, - mAj/i(mp) + m|/3(mp)], 
G 

^ycbVuDaifp[{mA, + mAj^i(mp) - mlgsiml)], 



(58) 



and 



T{A, ^ A,y) 
a(A, ^ A,\/) 



2(|5|^+|P2|') + ^(|5 + D|^ + |Pi|^) 



|P^|£;A.+mA, 
Stt mAi, [ 

AmlRe{S*P2) + 2E^Re{S + DfPi 
2mf.(|5|2 + |P2|2) + El{\S + DP + |Pi|2) ^ 



(59) 



where 



-A, 



IP' 



712 



Ev{Ea^ + tua. 



■{Ai - mA.As), 
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with 



\P'\ ,m\, + m\ 



^1 = --^VcbVuj:)aifvmv[gi{rr?y) + g2{ml,){m^^ - mAj], 



^2 = -2--/|V;fcVc;^ai/ymy5(2(m^), 



5i = -^VcbVuDaifvmvlfiirriy) - /2(my)(mA, + mAj], 
G 

^2 = 2--^V;bVf;^ai/vmy/2(my). (61) 



Then from Eqs.(p6D-(|6T]), we obtain the numerical resuhs for the decay widths 
and asymmetry parameters. In Table 2 we list the results for mr, = O.TOGeV. 
For other values of m^i, the results change only a little. The numbers without 
(with) brackets correspond to k = 0.02GeV^ (k = O.lOGeV'^). Again, the subscripts 
"0", "l/mg", and "l/mg + QCD" stand for the results without I/tuq and QCD 
corrections, with l/mq corrections, and with both l/rriQ and QCD corrections, 
respectively. In the calculations we have taken the following decay constants 

= 132MeV, fK = 156MeV, fn = 200MeV, fn, = 241MeV, 

/, = 216MeV, fK*=fp, fD = fD*, fD. = fDf 

Since the changes for the up-down asymmetries caused by l/rriQ and QCD cor- 
rections are very small, in Table 2 we only listed ai/mg+QCD- Furthermore, since to 
order ©(ttsA/mg) all the six form factors Fj, Gj {i = 1,2,3) can be expressed by 
one form factor, say Fi, which is canceled in a, the up-down asymmetries are model 
independent. Therefore, a does not depend on k. It can be seen from Table 2 that 
the predictions for the decay widths show a strong dependence on the parameters k 
in our model. In the future the experimental data will be used to fix this parameter 
and test our model. 
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Table 2: Predictions for the decay rates (in units 10^°s ^af, which is defined in 
Eq.(pl|)), and the asymmetry parameters for A;, AcP{V) 



Process 




ri/mg 


Tl/mQ+QCD 


"l/mq+QCD 


AO - 


A+TT- 


0.30 (0.56) 


0.36 (0.67) 


0.29 (0.55) 


-1.00 


Ar 


-A+p- 


0.44 (0.78) 


0.51 (0.94) 


0.42 (0.77) 


-0.89 




- A+D; 


1.03 (1.57) 


1.16 (1.81) 


1.02 (1.59) 


-0.98 


Ar 


- A+^r 


0.78 (1.17) 


0.89 (1.35) 


0.76 (1.15) 


-0.38 


Ar 


-y A+K~ 


0.022 (0.039) 


0.026 (0.048) 


0.021 (0.039) 


-1.00 


Ar 


-y A+K*- 


0.023 (0.041) 


0.027 (0.049) 


0.022 (0.040) 


-0.85 




A+D- 


0.037 (0.057) 


0.042 (0.066) 


0.036 (0.057) 


-0.98 


AO - 


A+D*- 


0.027 (0.041) 


0.031 (0.048) 


0.026 (0.040) 


-0.42 



In our previous workjl^, the Aj, Ac semileptonic and nonleptonic decay 



widths were calculated using a hadronic wave function model in the infinite momen- 
tum frame by combining the Drell-Yan type overlap integrals and the results from 
HQET to order 1 / mq. Comparing the results in our present B-S model with those in 
Refs. |ll3| , [T^ , we find that there is overlap between these two model predictions. The 
results with k = 0.02GeV'^ in the present model are close to those in Refs.[|l3], Q if 
the average transverse momentum of the heavy quark is chosen as 400MeV. 

The Cabibbo-allowed nonleptonic decay widths have also been calculated in the 
nonrelativistic quark model approach p9[, where the form factors are calculated at 



the zero-recoil point and then extrapolated to other u values under the assumption 
of a dipole behavior. It seems that the predictions in this model are close to those 
in our present work if we choose k = 0.02GeV^. 

V. Summary and discussion 

In the present work, we assume that a heavy baryon Ag is composed of a heavy 
quark, Q, and a scalar light diquark. Based on this picture, we analyze the l/mq 
corrections to the B-S equation for Aq which was established in the limit mg —>■ oo 
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in previous work ||Tl|. We find that in addition to the one B-S scalar function 
when niQ — » oo, two more scalar functions, (pip{p) and 4>2p{p), are needed at order 
l/rriQ. (p2p{p) is related to 4>op{p) directly [Eq.(|T^)j. Furthermore, with the aid 
of the reasonable constraint on the B-S kernel at order l/mq, Eg . (|20|) , and Luke's 
theorem, 0ip(p) can also be related to the B-S scalar function in the leading order. 
Hence we do not need to solve explicitly for 0ip(p) and (f)2p{p) any more. The B-S 
wave function in the leading order of I/ttlq expansion was obtained numerically 
by assuming the kernel for the B-S equation in the limit tuq — oo to consist of 
a scalar confinement term and a one-gluon-exchange term. On the other hand, 
all the six form factors for A;, Ac are related to each other to order l/mq, 
as indicated from HQET. We determine these form factors by expressing them in 
terms of the B-S wave functions. We also show explicitly that the results from our 
model are consistent with HQET to order l/mq. We also discuss the dependence 
of our numerical results on the various parameters in our model. It is found that 
Fj, Gi {i = 1,2,3) are insensitive to the binding energy, at order l/niq, and their 
dependence on the diquark mass, rrip), is mild. However, the numerical solutions are 
very sensitive to the parameter k. 

Furthermore, we apply our solutions for the weak decay form factors to calculate 
the differential and total decay widths for the semileptonic decays A^ AJu, and 
the nonleptonic decay widths for Af, — > AcP{V). The QCD corrections are also 
included, and found to be comparable with the l/rriq corrections. Again the numer- 
ical results for the decay widths mostly depend on k. We also compare our results 
with other models, including the hadronic wave function model and the norelativis- 
tic quark model, where l/mq corrections are also included. Generally predictions 
from these models are consistent with each other if we take into account the range of 
model parameters. Data from the future experiments will help to fix the parameters 
and allow one to test these models. 
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Besides the uncertainties from the parameters in our model, higher order correc- 
tions such as 0(1/771^) and 0{asA/mQ) will modify our results. However, we expect 
them to be small. Furthermore, we take a phenomenologically inspired form for the 
kernel of the B-S equation and use the covariant instantaneous approximation while 
solving the B-S equation. All these ansatze should be tested by the forthcoming 
experiments. 
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